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Abstract 

The thermodynamics of a long self-repulsive chain is studied. In D < 4 

dimensions it shows two distinctly different regimes, corresponding to weak 

and strong correlations in the system. A simple microscopic analytical model 

is presented which successfully describes both the regimes. The self-consistent 

scheme is used, in which the center of mass of a chain is fixed explicitly. This 

allows to take correlations into account in an indirect manner. 
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Problem of a crossover between the strong-correlation regime and weak correlations is 
of a general interest for the solid-state physics. A textbook example is here the theory of 
phase transitions Not too close to the transition point Landau theory is valid, yielding 
the dimension-independent Landau set of critical indices. On the other hand, correlations in 
a critical region in a very vicinity of the transition result in different (dimension-dependent) 
values of critical indices, as it can be described by the renormalisation-group analysis 
While the two regimes are well understood in itself, no general microscopic analytical formula 
describes them simultaneously. 

In this paper we draw a formula of this kind, which describes the behaviour of a simpler 
system, namely the self-repulsive chain. The peculiarity of this system is that it does not 
show a phase transition, but the two regimes mentioned above are still present. We believe 
that an absence of a spontaneous symmetry break makes things simpler, while the physics 
of a crossover remains somehow similar to more complicated problems like the description 
of the phase transition; for instance the problem of a self-interacting chain is closely related 
to the theory of percolation 

Let us first define the formal (mathematical) problem. We consider paths of a random 
walker on a cubic lattice in D dimensions, which returns to the starting point after L steps 
(see Figure 1 for the illustration); we suppose that L ^ 1. These paths are closed polygons 
(maybe, self-crossing); positions of their vertexes x-i fulfill the condition 

\\xi - Xi+i\ \ = 1; xo = XL- (1) 

We define the number of self-crosses for a given path X = {xo,xi, ...,xl-i} by the formula 

iV(X) = i5:5(x,;x,), (2) 

where 5 is a Kronecker symbol, and introduce a statistical average as follows: 

</> =^-^E/We"''''^''^; (3) 

X 

Z = Y,e-^N(,x)_ 

X 
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Here sums are taken over all paths X, / is a function of X, and /3 is a small positive 
parameter. In particular, we shall monitor the dispersion < {x — u)^ >, where u = J^i^i 
is the position of the center of mass of a chain. 

Physically, the above equations describe the statistics of a loop-chain with a potential 
energy equal to the number of its self-crosses N, at the inverse temperature f3. 

Since /5 -C 1, for a sufficiently small L the statistical weight e~^'^ is always close to 1, i.e. 
correlations in the system are negligible. Therefore we deal with a standard random-walk 
problem; for this case < {x — uY >oc L. We shall refer to this case as an uncorrelated 
regime. 

Let us establish a parameter which determines the magnitude of fluctuations. Since the 
typical volume occupied by the chain is proportional to L^/^ in the uncorrelated regime, 
the average number of its self-crosses can be estimated as X oc l?^^!"^. Therefore the 
dimensionless length A = 0^l^'^~^^L can be introduced; the uncorrelated regime occurs at 
/3X <^ 1, that corresponds to A ^ 1. 

In the opposite case of a large A, correlations in the chain affect its statistics strongly. 
The problem is similar to the well-studied case of the statistics of a self-avoiding polymer 
molecule [^^. The scahng laws are dimension- dependent and can be calculated by the 
renormalisation-group technique. For this "fractal" limit < [x — xx L!^ with = 2 in 
ID and takes non-integer values in 2D and 3D, which are respectively close to 4/3 and 6/5. 
In D > 4 the fractal regime is virtually absent, as takes the "uncorrelated" value 1 (there 
are logarithmic corrections in D = 4). 

As it was declared above, the present paper is aimed to establish a formula that describes 
(at least qualitatively) both the uncorrelated regime and the fractal regime. The outline is 
as follows. First, the continuous analog (|^ of the chain under study is introduced. Then, 
we switch to a kind of the self-consistent potential approximation (pi]). It yields simple 
formulae (p!5| - |T9D , that show the desired asymptotes for < [x — uf' >. The result is then 
compared with numerical calculations. 

Let us construct a continuous analog. One can neglect the short-range correlations, 

3 



as the critical behaviour in D < 4 is due to the long-range ones 0. Let r S> 1 be a 

"physically infinitesimal" length, on which correlations can be neglected; accordingly to 

the above estimations it should be r <^ Since P is small, in D < 4 these two 

inequahties can be fulfilled simultaneously. Consider the paths passing through some Xi. 

Then the distribution function of a; at i + r is given by 

, , / D(x — Xi)'^\ 
Pi+r{x; Xi) = po exp I — I , (4) 

where Po{t) is a normalisation factor. The number of paths passing simultaneously through 
given Xi and Xi+r, is obviously proportional to Pi+r{xi+r'i Xi). 

Introduce also the "partial" average number of self-crosses riij, which counts only the 
crosses between the "sub-chains" of the length r located near i and j. It is defined as 

TT-ij^Xif Xj-j-^ XiY ■3^{j+l)r; ■^jr '^{i + l)T) ^ ^ y '^(■^A.'i Xjyi) (5) 

fe,m 

where indices k and m run only within the intervals ... z(r + 1) — 1 and jr ... j(r + 
1) — 1, respectively, whereas the averaging is carried out over the paths passing through 

Xiri •^(i+l)T) Xjf, and XQ-i-I^t-. 

Full path L can be divided into L/t parts of the length r, so that, for example, partition 
function Z takes the form 

Z = Zo(r) E exp f- ^"^-7^'+^^-"' - f E . (6) 

Here the pre-exponential sum is taken over all possible sets {xq, Xt, xl-t}- 

Obviously, n{r,pi,p2) peaks near r = and falls fast at r ^ 00. It can be verified |^ 
that Y^r'^ij') Pi) P2) = at any pi,p2- So the continuous limit of (H) is as follows: 



Z = ZoJ [Dx] exp 



^^'^^^Idl - ^ S{x{l) - x{l'))dldl' ) . (7) 
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Here 6 is the Dirac 5-function. 

Self-parallel shifts of the whole chain are irrelevant for the statistical properties under 
study, so we can take the path-integral only over the closed paths with the center of mass 
fixed at zero: 



u 



[^xdl = 0; x{0)=x{L). (8) 

^0 



Continuous analog of the formula @ for statistical averages can be written in the same 
manner. In particular, the distribution function of x is equal to 

p{xo) = Zp{xo)/Z, (9) 

where Zp{xo) is a contribution to Z from the trajectories passing though the given point xq 
at / = 0, that is x(0) = Xq. Function p{xq) is normalised: 

p{xo)dxo = 1. (10) 

The key assumption of the present model is to approximate the second term in the 
exponent by the interaction with harmonic potential: 



-^11 S{x{l) - x{l'))dldl' ^ [^'+ 2 j^^' ^^^^ 
where Vt is defined in a self-consistent way: 

= -l3LD-\V^p%,=o (12) 

(the value of Uq may remain undefined, as it drops out of the final formulae). 

After the approximation (|1T|) is done, the distribution function can be found straight- 
forwardly. The problem resembles Feynman's path-integral expressions for the statistics of 
a quantum particle [0], with however an additional condition for trajectories / xdl = 0. Let 
us introduce a notation 

^-(0 = / i^-Uo- dl (13) 



/o V 2 2 ^ 

and define the "minimal- action" trajectory x{l; xq), which satisfy (H) and delivers a minimum 
of S at given x(0) = Xq. The conditional minimisation results in the problem of a classical 
oscillator in the external field 

|| + fi^x-/ = 0, (14) 



where the value of / should be chosen to fulfil the condition The following minimal 
value of S is obtained: 

Srmnix,) = -UoL + ^^^^^, (15) 

A(^) = — 



3 \ 1 — cuctgcij J 
uj = nL/2. 

Consider now possible deviations from this trajectory: x{l) = x{l)+x{l), where x(0) = and 
the conditions (H) are satisfied. As the potential is harmonic it appears that Sx — Sn2in~^Sxi 
and Zp can be factorised: 

Zp{xo) = Zexp{-Smin{xo)). (16) 



Here Z does not depend on xq. The last two formulae and expressions (P, |TOD allow to find 
out the distribution function: 

VM=[-^] exp( ^) (17) 



The dispersion is equal to 

< >= — (18) 
Substitution of V^p in formula ( |I2D gives the equation 

The left-hand and right-hand sides of this equation are schematically plotted in Figure 
2 as functions of uj. There is a single root within the range (0,7r). Its position is governed 
by A = I3'^/^*~'^^L. At A < 1 we obtain ^ and therefore K{uj) ^ 1, so 

<x^> Ia«i = ^. (20) 
The " uncorrelated" dependence < >oc L is reproduced; moreover the factor 1/12 is also 



correct. At the opposite limit ^^(A ^ 1) ^ vr, and it follows from ( p!q - |T9|) that 
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In particular < x"^ >oc L in 4D and < >oc in ID. The indices in 2D and 3D are also 
quite reasonable. In D > 4 the scheme is hardly suitable, as it is impossible to fulfill the 
inequality 1 <^ r <^ ^-2/(4-1)) g^j^^ ^j-^g continuum analog (|^) cannot be used. 

It is important to discuss the following peculiarity of the approximation made. Usually, 
the use of mean-field schemes means that correlations like < x(/i)x(/2) > at /i 7^ I2 are 
totally neglected. In our scheme however some correlations are still present, because the 
center of mass of the trajectory is fixed (formula (H)). Namely that circumstance results in 
a non-trivial dependence of Smm{L) and the dimension-dependent critical index in the last 
formula. 

We have also compared the result obtained with Monte-Carlo calculations in ID and 2D. 
Figure 3 shows bi-logarithmic plots for the dispersion < > versus L aX (3 = 10~^ in ID 
and (3 = 10~^ in 2D. Thin lines are asymptotes ( |20| , pTD . An agreement is quite good even 
quantitatively. 

The final remark is that although we have considered only a particular random-walk 
problem (|T]) with unitary steps, it is clear that the model can be applied for any random- 
walk self-repulsing chains. It is only required that (a) the repulsion is short-range and (b) 
the central limit theorem can be applied, yielding the Gaussian distribution @). Once these 
requirements are fulfilled, the continuous analog of the universal form (0) can be established. 

The author is grateful to P.A. Prudkovskii, A. A. Nikulin, and A. A. Fedyanin, who made 
valuable stylistic comments. The work was partially supported by "Scientific Schools" pro- 
gram (grant 96-1596476). 
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FIGURES 

FIG. 1. Examples of the self-crossing closed chains in D = 2. In the central chain there is a 

single self-cross (positions of vertexes 1 and 3 are the same). In the right chain there are 3 whose 
pairs (3-5, 3-7, and 5-7). 

FIG. 2. The sketch of A(f2) (thick line) and the right-hand side of the equation (19) (thin lines). 
Note that in D < 4 the later is always a monotonic increasing function of w, so there is always a 
single root of (19). 

FIG. 3. Monte-Carlo calculations of the dispersion < > in ID and 2D (points), compared 
with estimations from formula (18). Thin lines show the displacive and fractal asymptotes. 
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